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ABSTRACT: In this paper, we introduce an approach to increase density of field-effect transistors framework 
Schmitt trigger delay element. Framework the approach we consider manufacturing of the inverter in 
heterostructure with specific configuration. Several required areas of the heterostructure should be doped by 
diffusion or ion implantation. After that, dopant and radiation defects should by annealed framework by 
optimized scheme. We also consider an approach to decrease value of mismatch-induced stress in the 
considered heterostructure. We introduce an analytical approach to analyze the mass and heat transport in 
heterostructures during manufacturing of integrated circuits with account mismatch-induced stress. 
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INTRODUCTION 

In the present time several problems of the solid 
state electronics (such as increasing of performance, 
reliability and density of elements of integrated 
circuits: diodes, field-effect and bipolar transistors) 
are intensively solving by mathematical modeling [1- 
6]. To increase the performance of these devices, it is 
attracted an interest of the materials with higher 
values of charge carriers mobility [7-10]. One way to 
decrease the dimensions of elements of integrated 
circuits by manufacturing them in thin film 
heterostructures [3-5,11]. In this case, it is possible 
to use inhomogeneity of heterostructure and 
necessary to optimize and doping of electronic 
materials [12] and to develop epitaxial technology to 
improve these materials (including analysis of 
mismatch induced stress) [13-15]. An alternative 
approach to increase the dimensions of integrated 
circuits are using of laser and microwave types of 
annealing [16-18]. 


In the present study, we introduce an approach to 
manufacture field-effect transistors. This approach 
gives a possibility to decrease their dimensions with 
increasing their density of the frame work using 
Schmitt trigger delay element. We also consider 
possibility to decrease mismatch-induced stress to 
decrease the quantity of defects, generated due to 
stress. In this paper, we consider a heterostructure, 
which consist of a substrate and an epitaxial layer 
(see Fig. 1). We also consider a buffer layer between 
the substrate and the epitaxial layer. The epitaxial 
layer includes into itself several sections, which were 
manufactured by using other materials. These 
sections have been doped by diffusion or ion 
implantation to manufacture the required types of 
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conductivity (p or n). These areas became sources, 
drains and gates (see Fig. 1). After this doping, it is 
required annealing of dopant and/or radiation 
defects. The main aim of the article is to analysis the 
redistribution of dopant and radiation defects to 
determine conditions, which correspond to 
decreasing of elements of the considered delay 
element and at the same time to increase their 
density. At the same time we consider also a 
possibility to decrease mismatch-induced stress. 


Fig. 1a. Heterostructure with a substrate, epitaxial 
layers and buffer layer (view from top) and 
structure of the considered Schmitt trigger delay 
element [24] 
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To solve our ideas, we determine and analyzed 
spatio-temporal distribution of concentration of 
dopant in the considered heterostructure. We 
determine the distribution by solving the second 
Fick's law in the following form [1,19-23] 


_-Eptaxil Layer 


Fig. 1b. Heterostructure with a substrate, epitaxial 
layers and buffer layer (view from side) 
Method of solution 
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Here C(x,y,z,t) is the spatio-temporal distribution of concentration of dopant; Q is the atomic volume of dopant; Vs 
L, 

is the symbol of surficial gradient; (C(x, y,Z,t)d Z is the surficial concentration of dopant on interface 
0 


between layers of heterostructure (in this situation we assume, that Z-axis is perpendicular to interface between 
layers of heterostructure); 4u(x,y,z,t) is the chemical potential due to the presence of mismatch-induced stress; D 
and Ds are the coefficients of volumetric and surficial diffusions. Values of dopant diffusions coefficients depends 
on properties of materials of heterostructure, speed of heating and cooling of materials during annealing and 
spatio-temporal distribution of concentration of dopant. Dependences of dopant diffusions coefficients on 
parameters could be approximated by the following relations [21-23] 


C’ (x, y, Zt) V (x, y, Zt) V(x, y,Z,t) 
D, = D,\x, y, z,T )| 1+ € —>—— } 1 
G Be: Yk ) +5 b(n y,27) +O} V a?) vy 


C’ (x, y, z,t) V (x, y, z,t) V(x, y, z,t) 
1 td 
+85 P’ (x, y,z,T) + Cy ve +G5 (aa 


Here Di (x,y,z,T) and Dis (x,y,z,T) are the spatial (due to accounting all layers of heterostructure) and temperature 
(due to Arrhenius law) dependences of dopant diffusion coefficients; T is the temperature of annealing; P(x,y,z,T) 
is the limit of solubility of dopant; parameter y depends on properties of materials and could be integer in the 
following interval y €[1,3] [21]; V (xy,z,t) is the spatio-temporal distribution of concentration of radiation 
vacancies; V’ is the equilibrium distribution of vacancies. Concentration dependence of dopant diffusion 
coefficient has been described in details in [21]. Spatio-temporal distributions of concentration of point radiation 
defects have been determined by solving the following system of equations [19,22,23] 


D, = Dy, (x, y,z,T) (2) 
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I (xy,z,0)=ft (x,y,z), V (x,y,z,0)=fv (x,y,z). 
Here / (x,y,z,t) is the spatio-temporal distribution of concentration of radiation interstitials; I* is the equilibrium 
distribution of interstitials; Di(x,y,z,T), Dv(x,y,z,T), Dis(x,y, Z,T), Dvs(x,y,z,T) are the coefficients of volumetric and 
surficial diffusions of interstitials and vacancies, respectively; terms V2(x,y,z,t) and /2(x,y,z,t) correspond to 
generation of divacancies and diinterstitials, respectively (see, for example, [23] and appropriate references in 
this book); kiv(x,y,Z,T), kii(x,y,z,T) and kv,v(x,y,z,T) are the parameters of recombination of point radiation defects 
and generation of their complexes. 
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Spatio-temporal distributions of divacancies @ (x,y,z,t) and diinterstitials ® (x,y,z, t) could be determined by 
solving the following system of equations [19,22,23] 
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Here Dai(x,y,Z,T), Dov(X,y,2Z,T), Dos (%y,Z,T) and Dovs(x,y,z,T) are the coefficients of volumetric and surficial 
diffusions of complexes of radiation defects; ki(x,y,z,T) and kv(x,y,z,T) are the parameters of decay of complexes of 
radiation defects. 


Chemical potential su in Eq.(1) could be determine by the following relation [19] 
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M=E(Z)Qoy [Ui(Xy,2,0)+ui(Xy,2,)]/2, (7) 
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ui, uj are the components u:,(x,y,Z,t), uy(x,y,z,t) and u-(x,y,z,t) of the displacement vector u (se Y,Z, t ); Xi, Xjare the 


coordinate x, y, z. The Eq. (3) could be transform to the following form 


where E(z) is the Young modulus, oj is the stress tensor; Uj; = is the deformation tensor; 
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where ois Poisson coefficient; & = (as-az.)/az. is the mismatch parameter; ds, az are lattice distances of the 
substrate and the epitaxial layer; K is the modulus of uniform compression; / is the coefficient of thermal 
expansion; T; is the equilibrium temperature, which coincide (for our case) with room temperature. Components 
of displacement vector could be obtained by solution of the following equations [25] 
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Conditions for the system of Eq. (8) could be written in the form 
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We determine spatio-temporal distributions of concentrations of dopant and radiati-on defects by solving the 
Eqs.(1), (3) and (5) framework standard method of averaging of function corrections [26]. Previously we 
transform the Eqs.(1), (3) and (5) to the following form with account initial distributions of the considered 
concentrations 
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Farther we replace concentrations of dopant and radiation defects in right sides of Eqs. (1a), (3a) and (5a), and 
their not yet known average values ai, In this situation we obtain equations for the first-order approximations of 
the required concentrations in the following form 
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Integration of the left and right sides of the Eqs. (1b), (3b) and (5b) on time gives us possibility to obtain relations 
for above approximation in the final form 
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We determine average values of the first-order approximations of concentrations of dopant and radiation defects 
by the following standard relation [26] 
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Substitution of the relations (1c), (3c) and (5c) into relation (9) gives us possibility to obtain required average 
values in the following form 
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We determine approximations of the second and higher orders of concentrations of dopant and radiation defects 
framework standard iterative procedure of method of averaging of function corrections [26]. Framework this 
procedure to determine approximations of the n-th order of concentrations of dopant and radiation defects we 
replace the required concentrations in the Eqs. (1c), (3c), (5c) on the following sum anptp n1 (X,y,z,t). The 
replacement leads to the following transformation of the appropriate equations 
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Integration of the left and the right sides of Eqs. (1d), (3d) and (5d) gives us possibility to obtain relations for the 
required concentrations in the final form 
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Average values of the second-order approximations of required approximations by using the following standard 
relation [26] 
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Further, we have determined the solutions of Eqs.(8), ie. components of displacement vector. To determine the 
first-order approximations of the considered components framework method of averaging of function 
corrections we replace the required functions in the right sides of the equations by their not yet known average 
values a. The substitution leads to the following result 
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Integration of the left and the right sides of the above relations on time t leads to the following result 
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Approximations of the second and higher orders of components of displacement vector could be determined by 
using standard replacement of the required components on the following sums ait+ui({x,y,z,t) [26]. The 
replacement leads to the following result 
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Integration of the left and right sides of the above relations on time t leads to the following result 
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In this framework of the study, we determined the 
concentration of dopant, concentrations of radiation 
defects and components of displacement vector by 
using the second-order approximation framework 
method of averaging of function corrections. This 
approximation is usually enough good approximation 
to make qualitative analysis and to obtain some 
quantitative results. All obtained results have been 
checked by comparison with results of numerical 
simulations. 


DISCUSSION 

In this section, we analyzed dynamics _ of 
redistributions of dopant and radiation defects 
during annealing and under influence of mismatch- 
induced _ stress. _—‘ Typical distributions of 
concentrations of dopant in heterostructures are 
presented on Figs. 2 and 3 for diffusion and ion types 
of doping, respectively. These distributions have 
been calculated for the case, when value of dopant 
diffusion coefficient in the epitaxial layer is larger, 
than in the substrate. The figures show, that 
inhomogeneity of heterostructure gives us possibility 
to increase sharpness of p-n- junctions. At the same 


time one can find increasing homogeneity of dopant 
distribution in doped part of epitaxial layer. 
Increasing of sharpness of p-n-junction gives us 
possibility to decrease switching time. The second 
effect leads to decreasing local heating of materials 
during functioning of p-n-junction or decreasing of 
dimensions of the p-n-junction for fixed maximal 
value of local overheat. However framework this 
approach of manufacturing of bipolar transistor it is 
necessary to optimize annealing of dopant and/or 
radiation defects. Reason of this optimization is 
following. If annealing time is small, the dopant did 
not achieve any interfaces between materials of 
heterostructure. In this situation one cannot find any 
modifications of distribution of concentration of 
dopant. If annealing time is large, distribution of 
concentration of dopant is too homogenous. We 
optimize annealing time framework recently 
introduces approach [15,25-32]. Framework this 
criterion we approximate real distribution of 
concentration of dopant by step-wise function (see 
Figs. 4 and 5). Farther we determine optimal values 
of annealing time by minimization of the following 
mean-squared error. 
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Fig.2. Distributions of concentration of infused 
dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial 
layer substrate. Increasing of number of curve 
corresponds to increasing of difference between 
values of dopant diffusion coefficient in layers of 
heterostructure under condition, when value of 
dopant diffusion coefficient in epitaxial layer is 
larger, than value of dopant diffusion coefficient in 
substrate 
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Fig.3. Distributions of concentration of implanted 
dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial 
layer substrate. Curves 1 and 3 corresponds to 
annealing time © = 0.0048(L,?+Ly?+Lz*)/Do. Curves 2 
and 4 corresponds to annealing time © = 
0.0057 (Lx?+Ly*+Lz*)/Do. Curves 1 and 2 correspond to 
homogenous sample. Curves 3 and 4 corresponds to 
heterostructure under condition, when value of 
dopant diffusion coefficient in epitaxial layer is 
larger, than value of dopant diffusion coefficient in 
substrate 
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Fig.4 Spatial distributions of dopant in 
heterostructure after dopant infusion. Curve 1 is 
idealized distribution of dopant. Curves 2-4 are real 
distributions of dopant for different values of 
annealing time. Increasing of number of curve 
corresponds to increasing of annealing time 
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Fig 5. Spatial distributions of dopant in 


heterostructure after ion implantation. Curve 1 is 
idealized distribution of dopant. Curves 2-4 are real 
distributions of dopant for different values of 
annealing time. Increasing of number of curve 
corresponds to increasing of annealing time. 
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Fig.6. Dependences of dimensionless optimal 
annealing time for doping by diffusion, which have 
been obtained by minimization of mean-squared 
error, on several parameters. Curve 1 is the 
dependence of dimensionless optimal annealing time 
on the relation a/L and € = y= 0 for equal to each 
other values of dopant diffusion coefficient in all 
parts of heterostructure. Curve 2 is the dependence 
of dimensionless optimal annealing time on value of 
parameter ¢for a/L=1/2 and €= y= 0. Curve 3 is the 
dependence of dimensionless optimal annealing time 
on value of parameter & for a/L=1/2 and ¢ = y= 0. 
Curve 4 is the dependence of dimensionless optimal 
annealing time on value of parameter y for a/L=1/2 
and e= €=0 
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where yw (x,y,z) is the approximation function. 
Dependences of optimal values of annealing time on 
parameters are presented on Figs. 6 and 7 for 
diffusion and ion types of doping, respectively. It 
should be noted, that it is necessary to anneal 
radiation defects after ion implantation. One could 
find spreading of concentration of distribution 
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Fig.7.. Dependences of dimensionless optimal 
annealing time for doping by ion implantation, which 
have been obtained by minimization of mean- 
squared error, on several parameters. Curve 1 is the 
dependence of dimensionless optimal annealing time 
on the relation a/L and € = y= 0 for equal to each 
other values of dopant diffusion coefficient in all 
parts of heterostructure. Curve 2 is the dependence 
of dimensionless optimal annealing time on value of 
parameter ¢for a/L=1/2 and = y= 0. Curve 3 is the 
dependence of dimensionless optimal annealing time 
on value of parameter & for a/L=1/2 and ¢ = y= 0. 
Curve 4 is the dependence of dimensionless optimal 
annealing time on value of parameter y for a/L=1/2 
and = €=0 
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Fig. 8. Normalized dependences of component uz of 
displacement vector on coordinate z for nonporous 
(curve 1) and porous (curve 2) epitaxial layers 


of dopant during this annealing. In the ideal case 
distribution of dopant achieves appropriate 
interfaces between materials of heterostructure 
during annealing of radiation defects. If dopant did 
not achieve any interfaces during annealing of 
radiation defects, it is practicably to additionally 
anneal the dopant. In this situation optimal value of 
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additional annealing time of implanted dopant is 
smaller, than annealing time of infused dopant. 
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Fig. 9. Normalized distributions of charge carrier 
mobility in the considered heterostructure. Curve 1 
corresponds to the heterostructure, which has been 
considered in Fig. 1. Curve 2 correspond to a 
homogenous material with averaged parameters of 
heterostructure from Fig. 1 


Further we analyzed influence of relaxation of 
mechanical stress on distribution of dopant in doped 
areas of heterostructure. Under following condition 
&< 0 one can find compression of distribution of 
concentration of dopant near interface between 
materials of heterostructure. Contrary (at «>0) one 
can find spreading of distribution of concentration of 
dopant in this area. This changing of distribution of 
concentration of dopant could be at least partially 
compensated by using laser annealing [29]. This type 
of annealing gives us possibility to accelerate 
diffusion of dopant and another process in annealed 
area due to inhomogenous distribution of 
temperature and Arrhenius law. Accounting 
relaxation of mismatch-induced stress _ in 
heterostructure could leads to changing of optimal 
values of annealing time. Mismatch-induced stress 
could be used to increase density of elements of 
integrated circuits. On the other hand could leads to 
generation dislocations of the discrepancy. Fig. 8 
shows distributions of component of displacement 
vector, which is perpendicular to interface between 
layers of heterostructure. 


CONCLUSION 


In this paper, we model the redistribution of 
infused and implanted dopants with account 
relaxation mismatch-induced stress during 
manufacturing field-effect heterotransistors, which 
included in Schmitt trigger delay element. We 
formulate the recommendations for optimization of 
annealing to decrease dimensions of transistors and 
to increase their density. We suggested the 
recommendations to decrease mismatch-induced 
stress. Analytical approach to model diffusion and ion 
types of doping with account concurrent changing of 
parameters in space and time has been introduced. 
At the same time this approach gives us possibility to 
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take into account nonlinearity of considered 

processes. 
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